Differences in Coformer Interactions of the 2,4-Diaminopyrimidines Pyrimethamine and Trimethoprim

The identification and study of supramolecular synthons is a fundamental task in the design of pharmaceutical cocrystals. The malaria drug pyrimethamine (pyr) and the antibiotic trimethoprim (tmp) are both 2,4-diaminopyrimidine derivatives, providing the same C–NH2/N=C/C–NH2 and C–NH2/N=C interaction sites. In this article, we analyze and compare the synthons observed in the crystal structures of tmp and pyr cocrystals and molecular salts with sulfamethazine (smz), α-ketoglutaric acid (keto), oxalic acid (ox), sebacic acid (seb), and azeliac acid (az). We show that the same coformer interacts with different binding sites of the 2,4-diaminopyrimidine ring in the respective tmp and pyr cocrystals or binds at the same site but gives H bonding patterns with different graph set notions. Pyr·smz·CH3OH is the first crystal structure in which the interaction of the sulfa drug at the C–NH2/N=C/C–NH2 site with three parallel NH2···N, N···NHsulfonamide, and NH2···O=S H bonds is observed. The main synthon in (tmp+)(keto–).0.5H2O and (tmp+)2(ox2–)·2CH3OH is the motif of fused R21(6) and R12(5) rings instead of the R22(8) motif typically observed in tmp+ and pyr+ carboxylates. Tmp/az is a rare example of cocrystal-salt polymorphism where the two solid-state forms have the same composition, stoichiometry, and main synthon. Theoretical calculations were performed to understand the order of stability, which is tmp·az cocrystal > (tmp+)(az–) salt. Finally, two three-component tmp/sulfa drug/carboxylate cocrystals with a unique ternary synthon are described.


1
and R 1 2 (5) rings instead of the R 2 2 (8) motif typically observed in tmp + and pyr + carboxylates. Tmp/az is a rare example of cocrystal-salt polymorphism where the two solid-state forms have the same composition, stoichiometry, and main synthon. Theoretical calculations were performed to understand the order of stability, which is tmp·az cocrystal > (tmp + )(az − ) salt. Finally, two three-component tmp/sulfa drug/carboxylate cocrystals with a unique ternary synthon are described.

■ INTRODUCTION
The development of cocrystals is largely based on crystal engineering principles. In crystal engineering, supramolecular synthons between functional groups are used to manipulate the intermolecular interactions. By controlling the crystal packing, cocrystals with the desired properties can be designed. The largest and most widely studied class of cocrystals are pharmaceutical cocrystals, that is, two drugs or one drug and a pharmaceutically acceptable coformer held together in the same crystal lattice by noncovalent interactions. 1−4 Cocrystallization allows for the manipulation and optimization of the physicochemical properties of an active pharmaceutical ingredient (API) such as chemical and physical stability, hygroscopicity, solubility, and dissolution behavior without the need for chemically modifying the API molecule. APIs often have more than one functional group that can interact with coformers. Knowledge of the hierarchy of synthons is crucial to the rational design of new cocrystals. Furthermore, the presence of different functional groups or interaction sites can be exploited in the synthesis of higher-order cocrystals. 5−17 Cocrystals containing more than two components are significantly more difficult to isolate than binary cocrystals because not only the supramolecular synthons but also the size, shape, and solubility of the three components must match and there is always the risk that a more stable binary cocrystal crystallizes instead.
We have recently studied the formation of binary and ternary cocrystals of the malaria drug pyrimethamine (pyr, Figure 1). 18 The 2,4-diaminopyrimidine ring of pyr provides two distinct interaction sites for coformers, the C2−NH 2 /N3/ C4−NH 2 (donor−acceptor−donor, DAD) site and the C2− NH 2 /N1 (donor−acceptor, DA) site. We have investigated combinations of various ADA and AD coformers and carried out competition experiments with different AD coformers such as saccharin and monocarboxylic acids. We have now included dicarboxylic acid coformers and obtained more binary and ternary cocrystals of pyr with new structural motifs. The identification of new synthons and the understanding of synthon preferences is a fundamental objective of crystal engineering. We also report a cocrystallization study of the related 2,4-diaminopyrimidine derivative trimethoprim (tmp, Figure 1). Tmp is an antibiotic mainly used for the treatment of urinary tract infections. A few cocrystals and salts of tmp have been reported previously. 19−26 The crystal structure determination and the comparison of pyr-and tmp-coformer systems described in our present study revealed several interesting observations. (1) In contrast to pyr that formed a 1:1 salt with azeliac acid (az) as the stable form, the tmp/az system is one of the few examples of salt-cocrystal polymorphism. Polymorphism in multicomponent systems and the crystallization of solid-state forms with different ionization states are well established. However, systems with the same composition, stoichiometry, and the same main synthon that differ only in the location of the acidic proton at the same temperature are rare. 27−30 (2) While the sulfa drug sulfamethazine (smz) interacts with the DA binding site of tmpas observed for all previously reported tmp-sulfa drug cocrystals 19−23 a new synthon between smz and the DAD site of the 2,4-diaminopyrimidine ring is found in pyr·smz, namely, three parallel NH 2 ···N, N···NH sulfonamide , and NH 2 ··· OS H bonds. Oxalic acid (ox) and α-ketoglutaric acid (keto) also form different H bonding motifs with the 2,4diaminopyrimidine rings of tmp and pyr (fused R 2 1 (6) and R 1 2 (5) rings vs R 2 2 (8) rings). (3) Contrary to our expectation based on (2) that the cocrystallization with both a carboxylic acid (AD coformer) and a sulfa drug (ADA coformer) would give a ternary cocrystal with the two coformers occupying the two 2,4-diaminopyrimidine binding sites, crystal structure determination of (tmp + )(pim − )·stz and (tmp + )(seb − )·stz· 2H 2 O·C 3 H 6 O (pim = pimelic acid, seb = sebacic acid, stz = sulfathiazole) revealed a unique ternary synthon in which the SO and imido NH of the sulfa drug H bond to the amino nitrogen and carboxylate oxygen of the tmp + ···carboxylate pair. The reasons for the differences in the crystallization behavior of the closely related tmp and pyr molecules are discussed. ■ EXPERIMENTAL SECTION Materials. pyr, tmp, keto, az, pim, and seb were purchased from Tokyo Chemical Industry, Europe. ox and smz were obtained from Sigma-Aldrich. Sulfathiazole (stz) was obtained from Fluka Analytical. The solvents acetonitrile, methanol (Merck Millipore), ethyl acetate (Sigma-Aldrich), and acetone (Fisher Chemical) were of analytical grade and were used as received.
Solution Crystallization. 25 mg of pyr or tmp and 1 mol equiv of the respective coformer(s) were dissolved in the minimum amount of solvent, and the solvent was allowed to slowly evaporate from an open vial. Crystallization experiments were carried out in methanol, acetonitrile, acetone, and ethyl acetate. Within a few days, X-ray suitable single crystals of (pyr + ) 2 (ox 2− )·1.5H 2 O, (tmp + ) 2 (ox 2− )· 2CH 3 OH, (tmp + ) 2 (az 2− )·tmp·6H 2  Ball-Milling. Equimolar mixtures of tmp or pyr and the respective coformer(s) (120−150 mg in total, Table S1) were placed in 2 mL Eppendorf tubes. 20 μL of solvent and one 5 mm diameter stainless steel ball were added to each sample. The Eppendorf tubes were placed in an in-house 3D-printed six-tube sample holder. An oscillatory ball mill (Mixer Mill MM400, Retsch GmbH & Co., Germany) was used to mill the samples at 25 Hz for 20 min. The milled powder samples were analyzed immediately by X-ray powder diffraction.
Crystallization by Sublimation. The salt (pyr + )(az − ) form II was crystallized from the gas phase using an in-house sublimation apparatus. 31 Equimolar ratios of pyr and az were sublimed from both ends of a standard 15 × 160 mm test tube sealed under vacuum. Two heaters were used to sublime the components, and the temperatures at the pyr and az ends were set at 193.5 and 148.3°C, respectively. Pyr and az formed small crystals of (pyr + )(az − ) form II in the middle of the test tube after 17 h.
Slurry Experiments. Competitive slurry experiments were carried out for (tmp + ) 2 (az 2− )·tmp·6H 2 O and tmp·az to determine the most stable polymorph. 20 mg of each sample was mixed with 1 mL of methanol in a 10 mm diameter glass vial and the mixture was stirred for 48 h at room temperature, followed by drying in a vacuum oven at 40°C for 6 h and recording the powder X-ray pattern.
Differential Scanning Calorimetry. Differential scanning calorimetry (DSC) was carried out in open aluminum crucibles using an STA625 thermal analyzer (Rheometric Scientific, Piscataway, New Jersey). The DSC plots were recorded between 20 and 300°C with a heating rate of 10°C/min. Nitrogen was purged in ambient mode, and an indium standard was used for calibration.

■ RESULTS AND DISCUSSION
Tmp·pyr Cocrystal. The stable polymorphs of pyr and tmp have the same hydrogen-bonding motifs with pairwise H bonding at both amino-pyrimidine sites ( Figure 2a). 40,41 Pyr and tmp are therefore complementary coformers, and a cocrystal was indeed obtained from methanol. The crystal structure of the cocrystal of tmp and pyr is shown in Figure 2b. The asymmetric unit contains one tmp, one pyr, and one water molecule of crystallization. The R 2 2 (8) homosynthons with pairs of C2−NH 2 ···N1 and C4−NH 2 ···N3 H bonds in the Xray structures of pyr and tmp are replaced by the corresponding pyr···tmp heterosynthons. The second C4amino proton of tmp forms a bifurcated H bond with two methoxy groups of an adjacent tmp. In addition, one methoxy group of tmp acts as a H bond acceptor for the C2-amino group of pyr. Overall, pyr and tmp are linked to an infinite 2D structure. The water molecule of crystallization is disordered over three positions and interacts with the C4-amino group of pyr and with the C2-amino group of tmp.
The cocrystal can also be prepared by ball-milling an equimolar mixture of tmp and pyr in the presence of a trace amount of ethanol. The match of the XRPD pattern of the milled sample and the theoretical pattern calculated from the single crystal data of tmp·pyr·H 2 O is shown in Figure S1. The DSC plot of the cocrystal ( Figure S2) shows a dehydration endotherm at 124.1°C (weight loss in the TGA 3.1%; calcd. 3.2%) and a melting endotherm at 195.9°C. The latter temperature is lower than the melting points of tmp (200°C) and pyr (233°C).
Interaction of Dicarboxylic Acid Coformers with the DA Binding Site of Tmp and Pyr. A vast majority of crystal structures of cocrystals of tmp and pyr with a carboxylic acid coformer described in the literature show transfer of the carboxyl proton to the N1 nitrogen of pyr or tmp and the interaction of the carboxylate group with the N1−H + /C2− NH 2 site via a pair of parallel H bonds creating the R 2 2 (8) motif ( Figure 3). 43−49 In the following, we compare the H bonding interactions of the 2,4-diaminopyrimidine pharmacophore and carboxylic acid coformers for four carboxylic acids for which  Contrary to what one might expect, differences in the main H bonding motif were observed between tmp/keto and pyr/keto and between tmp/ox and pyr/ox, while tmp/az and pyr/az and tmp/seb and pyr/seb were found to differ in the ionization states of the coformers (Table 1). (Pyr + )(keto − ) 50 and (pyr + ) 2 (ox − )(ox 2− ) 0.5 ·2CH 3 OH 51 have the typical R 2 2 (8) motif; however, this motif is absent in the corresponding tmp salts. For the sake of simplicity, we use the notation pyr· X/tmp·X for a cocrystal and pyr + X − /tmp + X − for a salt with proton transfer from the coformer to the N1 nitrogen of pyr/ tmp.
The 2,4-diaminopyrimidine rings in pyr and tmp have almost the same pK a value (7.16 and 7.34, respectively), indicating very similar electronic properties and H bonding propensity. The fact that the R 1 2 (5)−R 2 1 (6) H bonding pattern is only observed in tmp + oxalates may be due to a steric hindrance effect of the ethyl substituent at C6 of pyr ( Figure  4a).
Milling of tmp and ox in a 2:1 ratio in the presence of traces of methanol gave an XRPD pattern that is a good match with the theoretical pattern of (tmp + ) 2 (ox 2− )·2CH 3 OH ( Figure S6). (Tmp + ) 2 (ox 2− )·6.5H 2 O and (pyr + ) 2 (ox 2− )·1.5H 2 O could not be obtained as phase-pure solids in bulk quantities, neither by    The tmp·az single crystal was obtained from acetonitrile. A comparison of the XRPD pattern of the crystalline sample isolated from acetonitrile with the theoretical powder pattern of tmp·az confirmed that the single crystal structure was representative of the bulk composition ( Figure S8). The single crystal of the (tmp + )(az − ) salt was isolated from ethyl acetate. However, XRPD analysis of the bulk sample showed the pattern of the cocrystal, indicating that the salt is a minor side product or a transient form. Attempts to prepare bulk quantities of the salt by rapid solvent evaporation were unsuccessful.
A third form, the ionic cocrystal (tmp + ) 2 (az 2− )·tmp·6H 2 O, crystallized from methanol (Figure 6d). The asymmetric unit comprises two tmp + cations (denoted A and B), one neutral tmp and one az 2− dianion and six water molecules of crystallization. In contrast to the 1:1 tmp·az cocrystal, (tmp + ) 2 (az 2− )·tmp·6H 2 O does not contain the C−NH 2 / N(H)···(H)OOC R 2 2 (8) motif. The az 2− dianion participates in H bonding with the C2-and C4-amino groups of tmp, the C2-amino group of tmp + A, and the C2-amino group of tmp + B and with the water molecule of crystallization. Tmp + A and tmp interact via a pair of H bonds between the N3/C4−NH 2 site of the cation and the C2−NH 2 /N3 site of the neutral molecule. There is also extensive H bonding between the water of crystallization and the amino groups of tmp + A, tmp + B, and tmp, a methoxy group of tmp, and N1 of tmp. Figure S9 shows the XRPD patterns after milling a 1:1 mixture of tmp and az in the presence of traces of ethanol, acetonitrile, and ethyl acetate. In all cases, the Bragg peaks of the cocrystal were observed. When a 3:1 mixture was milled in the presence of traces of H 2 O, the XRPD pattern of (tmp + ) 2 (az 2− )·tmp·6H 2 O was obtained ( Figure S10). On slurrying in methanol, (tmp + ) 2 (az 2− )·tmp·6H 2 O converted to tmp·az ( Figure S11). The DSC plot of (tmp + ) 2 (az 2− )·tmp· 6H 2 O shows endotherms at 90.4, 155.3, and 185.8°C ( Figure  S12a). The first endothermic event is accompanied by a 10.1% weight loss in the TA plot and can be assigned to the loss of the water of crystallization (calcd. 9.3%). The peak at 155.3°C corresponds to the melting point of tmp·az ( Figure S12b). When a sample of (tmp + ) 2 (az 2− )·tmp·6H 2 O was heated to 90°C under vacuum for 6 h, XRPD analysis showed the pattern of free tmp and minor peaks that may be assigned to the cocrystal ( Figure S13).
Stable Form of Tmp/Az and Pyr/Az. The only form that was identified in cocrystallization experiments for pyr and az is the salt (pyr + )(az − ). By contrast, solution crystallization of tmp/az suggests that the cocrystal is in the thermodynamically stable form, which was confirmed by slurry experiments. A slurry normally gives the thermodynamically stable product, and slurrying has been used in the literature to identify the most stable polymorph. 52,53 No change in the XRPD pattern was observed, when the cocrystal was slurried for 48 h in methanol ( Figure S14).
It has been reported for the ethionamide/salicylic acid system that proton transfer is not influenced by the nature of the solvent (polar/apolar; protic/aprotic) but that the crystallization of the salt/cocrystal polymorphs follows Ostwald's rule of stages. 27 Furthermore, we have previously shown that proton transfer can take place in prenucleation clusters in the absence of solvent, 54 which further corroborates the lack of a solvent effect on the ionization state of the coformers in the solid state. The intermediate formation of the (tmp + )(az − ) salt as a metastable transient form may be the reason why it was not possible to isolate bulk quantities.
The pK a values of the N1 nitrogen of pyr and for the deprotonation of the first carboxyl group of az are 7.34 and 4.15, respectively, so that salt formation for the pyr/az system is in line with the ΔpK a rule that states that a salt is usually obtained when the difference in the pK a of the two coformers is greater than 4, while a ΔpK a of <0 results in a cocrystal. 55,56 In the ΔpK a range 0 < pK a < 4, there is a salt-cocrystal continuum and a prediction of salt or cocrystal formation is not possible. The pK a value of tmp is 7.16 and is thus very close to that of pyr. Cruz-Cabeza has developed a model that estimates the probability for cocrystal and salt formation as 57 Using these equations, there is a 21% probability of tmp and az crystallizing as a cocrystal and a 79% probability for salt formation. The position of the proton is determined by the packing and the resulting lattice energy. While the amount of proton transfer can be affected by the crystalline environment, ΔpK a represents a useful predictive tool outside the continuum region. 56 Childs et al. compared the unintended application of pK a values to predict the ionization state in the solid state with the use of van der Waals radii, which were originally intended for the calculation of molecular volume but are now routinely used to estimate intermolecular distances in adjacent molecules. 56 Tmp·az and (tmp + )(az − ) are a rare example of cocrystal−salt polymorphism. The protonation state of the coformers in multicomponent crystals can vary for different stoichiometries or solvates. 58−60 However, only a small number of systems have been reported, where a salt or a cocrystal of the same composition and stoichiometry is obtained depending on the crystallization conditions; smz/ saccharine (smz/sac), 28,29 four amino acids/tartaric acid of which β-alanine/tartaric acid (bal/tar), has been examined in detail, 61,62 ethionamide/salicylic acid (eth/sal), 27 and a range of haloaniline/3,5-dinitrobenzoic acid systems. 30 The position of the proton may also be temperature-dependent. 63,64 Bal/ tar 62 and smz/sac 28,29 crystallize stochastically as salts or cocrystals from the same solution. Concomitant growth of salt and cocrystal polymorphs has also been observed for haloaniline/3,5-dinitrobenzoic acid, 30 and a metastable salt was identified as a transient form to the polymorphs of the 2:1 isonicotinamide−citric acid cocrystal. 65 The stable form, salt or cocrystal, was deduced in each case using crystal slurries and other methods, and in the case of the haloaniline/3,5dinitrobenzoic acid series, the energy differences per formula Crystal Growth & Design pubs.acs.org/crystal Article unit between the salt and cocrystal were calculated by periodic boundary DFT methods using the program VASP. In that series of compounds, the salts were more stable than the cocrystal forms without exception. In the case of 4-bromo-2methylaniline/3,5-dinitrobenzoic acid (brma/dnb), the salt was reported to be 8.5 kJ/mol more stable than the cocrystal. We have calculated the energy differences between the salt and the cocrystal for tmp/az and for the other systems described in the literature (Table 2), and in every case, periodic DFT calculations predict the stable form found by the experiment. Positive ΔE (sa-cc) values indicate that the cocrystal is more stable than the salt.
Crystal density is an indication of the efficiency of crystal packing and can be used as an indication of crystal stability. For example, crystal density has been used to select the most likely candidates in crystal structure prediction software. 66,67 In Table 2, crystal density correctly predicts the stable form in three out of the five examples, with bal/tar and brma/dnb having stable forms with lower densities than their metastable forms. Thus, at least in the examples in Table 2, crystal density is less effective than periodic DFT calculations in predicting the stability of crystal structures.
There is also one interesting difference between tmp/az and the other systems in Table 2 in that it was necessary to apply nonlinear distance restraints to the O−H groups H-bonded to nitrogen to prevent proton transfer. This of course does not imply that the crystallized salt is more stable than the cocrystal, rather this proton transfer suggests that there may be another salt form, which would crystallize in space group P-1 if suitable conditions for its isolation could be found.
Noteworthily, among the salt-cocrystal polymorph pairs listed in Table 2, tmp/az has the largest ΔpK a difference between the two coformers. To understand the higher stability of the tmp·az cocrystal and the difference to the pyr/az system, we have examined the packing of the tmp salt, tmp cocrystal, and pyr salt. The densities of (tmp + )(az − ), tmp·az, and (pyr + )(az − ) are 1.223, 1.307, and 1.252 g cm −3 , respectively, that is, the density of (pyr + )(az − ) is intermediate between those of tmp·az and (tmp + )(az − ). (Pyr + )(az − ) contains Cl and would be expected to have a higher density than either of the tmp/az forms. Using the Void program in the Oscail package, 69 the packing indices for tmp·az, (tmp + )(az − ), and (pyr + )(az − ) were calculated as 68.5, 70.5, and 63.7%. The Void program also finds a 3.3% void, which is large enough to insert one water molecule into the (pyr + )(az − ) structure. The program did not find any void in either of the tmp/az structures. The void locations in (pyr + )(az − ) are shown in the packing view down the c axis ( Figure S15). It is possible that the restriction imposed by the biphenyl linkage in the pyr molecule, which prevents the ring packing in a more efficient flat geometry, leads to a lower packing efficiency than is observed in the tmp/az structures, where the flexible nature of tmp may aid efficient packing. The observation of a salt rather than a cocrystal for pyr/az may be due to an extra lattice energy component provided by the ionic attractions in the salt. It is worthy to note that in tmp·az, both pyrimidine-N/NH 2 sites form a pair of H bonds with az-COOH, while in (tmp + )(az − ), only the N1/C2−NH 2 site forms the R 2 2 (8) motif and N3 is not involved in H bonding.
Different Interactions of Tmp and Pyr with Sulfa Drugs. The cocrystallization of tmp with the sulfa drugs sulfamethoxazole (smx), sulfametrole (smt), sulfamethoxypyridazine (smp), and smz is described in the literature, and the X-ray structures of the salts (tmp + )(smx − ), 19 (tmp + )(smt − ), 20 and (tmp + )(smp − )·1.5H 2 O 21 and of the cocrystals tmp·smp, 21 tmp·smz·CH 3 OH, 22 and tmp·2smz·H 2 O 23 have been reported. The sulfa drugs have an N(heterocycle)CH−NH-(sulfonamide) functionality, and in all cases, the heterocyclic nitrogen and the sulfonamide nitrogen form a pair of H bonds with the C2−NH 2 /N1 site of tmp (R 2 2 (8)). Interestingly, we obtained the single crystal structure of the smz cocrystal of pyr that showed binding of the sulfa drug via the DAD site of pyr ( Figure 7). In pyr·smz·CH 3 OH, one of the sulfonyl oxygens, the sulfonamide nitrogen and the ring nitrogen of smz form H bonds with the amino group at C4, the N3 nitrogen, and the amino group at C2 of pyr. The second sulfonyl oxygen accepts a H bond from the amino group of an adjacent smz. The methanol molecule of crystallization participates in H bonding with N1 of pyr and the amino group of smz. Again, the steric hindrance of the ethyl substituent on the carbon adjacent to N1 may be the reason why smz does not bind to the C2− NH 2 /N1 site as is the case in the tmp·smz cocrystal. It appears that the methanol of crystallization stabilizes the structure by filling the voids resulting from the formation of the ADA··· DAD synthon.
The DSC plot of pyr·smz ( Figure S16) shows an endotherm at 111.5°C followed by an exotherm at 136.2°C that can be assigned to the loss of the hydrogen-bonded methanol of recrystallization and the structural rearrangement of the desolvated form. The evaporation of methanol is confirmed  Ternary Cocrystals of Tmp. The observation that smz can act as an ADA coformer at the C2−NH 2 /N3/C4−NH 2 site of the 2,4-diaminopyrimidine ring in pyr·smz·CH 3 OH prompted us to explore the formation of ternary cocrystals of tmp/pyr, a sulfadrug, and a carboxylic acid. We expected that in such ternary mixtures, the carboxylic acid would have the highest affinity for the C2−NH 2 /N1 site of tmp/pyr: cocrystals of sulfa drugs and carboxylic acids are rather scarce. Very recently, it is has been stated in the literature that "the chance to form a successful cocrystal containing the sulfa drug and benzoic acid or one of its derivatives is low". 70 In the same article, a cocrystal screen of sulfathiazole (stz), sulfamethoxazole, and sulfapyridine with a range of coformers including carboxylic acids was reported. Out of the 30 sulfa drug−carboxylic acid combinations tested, only three were successful. These include the 3,5-dinitrobenzoate salt of stz, the sulfamethoxazole-3,5dinitrobenzoic acid cocrystal, and the 2-chloro-4-nitrobenzoate salt of sulfapyridine. A slightly higher success rate was reported for a smz cocrystal screen that led to the structural characterization of 7 smz cocrystals with aliphatic and aromatic carboxylic acids. In all cases, the carboxylic acid formed a pair of H bonds with the sulfonamide/N(heterocycle) site of either the amidine or the imidine tautomer of smz. 71 This R 2 2 (8) motif was also observed in the crystal structures of salts and cocrystals with benzoic acid, 72 salicylic acid, 73 anthranilic acid, 4-aminobenzoic acid, 74 4-aminosalicylic acid, acetylsalicylic acid, 75 2-nitrobenzoic acid, 76 4-nitrobenzoic acid, 77 2,4dinitrobenzoic acid, indole-2-carboxylic acid, 78 4-chlorobenzoic acid, 79 and 5-nitrosalicyclic acid. 80 By contrast, stz/3,5dinitrobenzoic acid, sulfamethoxazole/3,5-dinitrobenzoic acid, and sulfapyridine/2-chloro-4-nitrobenzoic acid all have different synthons. 70 In the stz·glutaric acid 81 and stz·4-nitrobenzoic acid 82 cocrystals, stz and the carboxylic acid both form homodimers. In stz·glutaric acid, the dimers are linked via a bifurcated H bond to the amino group and sulfonyl oxygen of two adjacent stz molecules. In stz·4-nitrobenzoic acid, the homodimers interact with each other via π−π contacts between the thiazole and nitrobenzene rings.
The carboxylic acid and sulfonamide combinations used in the screening study are listed in Tables S6 and S7. The isolation of a ternary cocrystal was only successful in the case of stz and tmp, and instead of binding to the DAD site of tmp, a new ternary synthon was observed with stz interacting with the tmp + ···carboxylate dimer via SO···NH 2 and NH···O H bonds. It is worthy to note that while in several cases, the binary pyr + /tmp + carboxylate salt crystallized from solution in the screening study for ternary cocrystals, no binary sulfa drug/ carboxylic acid cocrystal was obtained, which further confirms the low affinity of the carboxyl group for the sulfonamide/Nheterocycle binding site.
The H bonding patterns of the two ternary ionic cocrystals obtained, (tmp + )(pim − )·stz and (tmp + ) 2 (seb 2− )·stz·2H 2 O· C 3 H 6 O, are shown in Figure 8. In both structures, stz exists in its imido tautomeric form and forms the same ternary synthon with tmp + and the carboxylate group. H bonding between a sulfonyl oxygen and the amino group at C2 and between the proton on the heterocyclic ring nitrogen and a carboxylate oxygen connect the sulfa drug to the tmp + ··· carboxylate synthon, leading to a motif of fused R 3 2 (10) and R 2 2 (8) rings. Furthermore, both structures have the C4−NH 2 ··· N3 homosynthon. Additional hydrogen bonding of the C2and C4-amino groups with SO creates R 3 2 (8) rings. In (tmp + )(pim − )·stz, two pim monoanions dimerize into a 20membered ring, while the seb 2− dianion in (tmp + ) 2 (seb 2− )·stz· 2H 2 O·C 3 H 6 O lies on an inversion center, so that both carboxylate groups interact with stz and tmp + . It is interesting that in contrast to (tmp + )(pim − )·stz, no proton transfer from the carboxyl group to the 2,4-diaminopyrimidine nitrogen Crystal Growth & Design takes place when tmp and pim cocrystallize to give the binary cocrystal tmp·pim·0.5CH 3 CN ( Figure S17).
Milling experiments were carried out to investigate if the two ternary cocrystals could also be prepared mechanochemically. The XRPD patterns of the milled tmp/stz/seb sample (1:1:1 molar ratio) showed a new set of peaks and unreacted stz ( Figure S18). The new peaks are identical to the pattern obtained when a binary 1:1 mixture of tmp and seb is milled. In the case of the 1:1:1 tmp/stz/pim sample, the XRPD pattern matched the theoretical pattern of (tmp + )(pim − )·stz ( Figure S19), indicating that the ternary cocrystal can also be prepared mechanochemically.

■ CONCLUSIONS
Despite the identical H bonding functionality and almost identical pK a values, the 2,4-diaminopyrimidine groups in tmp and pyr do not necessarily form the same H bonding motif with a given coformer. Crystal structures have been obtained, where the same coformer binds to the N1/C2−NH 2 site of tmp (pairwise H bonding) but to the C2−NH 2 /N3/C4−NH 2 site of pyr (three parallel H bonds) or where the coformer interacts at the same site but gives H bonding patterns with different graph set notations, specifically fused R 2 1 (6) and R 1 2 (5) rings versus the R 2 2 (8) motif in carboxylates. These examples demonstrated the significant effect of the ethyl substituent on the C6 carbon of pyr. The tmp/az system shows salt-cocrystal polymorphism with the salt, presenting a transient metastable form, contrary to what one would expect on the basis of the ΔpK a rule, the Coulomb attraction contribution to the lattice energy, and the calculated packing indices of the salt and cocrystal forms. In ternary cocrystals, two coformers usually bind to two different functional groups or binding sites of the third component. Instead, all three components are connected in an R 3 2 (10) motif, generating a unique ternary motif in (tmp + )(pim − )·stz and (tmp + ) 2 (seb 2− )·stz·2H 2 O·C 3 H 6 O. We have reported three-component ionic cocrystals of pyr with ADA-and AD coformers in a previous paper. 18 The ADA-and AD-coformers bound to different pyr molecules in the crystal lattice. Apparently, the simultaneous interaction of two different coformers at both binding sites of the 2,4diaminopyrimidine ring is not favorable.

■ ACKNOWLEDGMENTS
This publication has emanated from research supported in part by a research grant from Science Foundation Ireland (SFI) and is co-funded under the European Regional Development Fund under grant number 12/RC/2275. The Irish Centre for High-End Computing (ICHEC) is thanked for the provision of computational resources under projects nuig02 and ngche102c.